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FULLBACK EXPONENTIAL ATTRACTOR FOR A 
CAHN-HILLIARD-NAVIER-STOKES SYSTEM IN 2D 

STEFANO BOSIA AND STEFANIA GATTI 

Abstract. We consider a model for the evolution of a mixture of two incompressible 
and partially immiscible Newtonian fluids in two dimensional bounded domain. More 
precisely, we address the well-known model H consisting of the Navier-Stokes equation 

^o , with non-autonomous external forcing term for the (average) fluid velocity, coupled with 

t . ' a convective Cahn-Hilliard equation with polynomial double-well potential describing 

f^ I the evolution of the relative density of atoms of one of the fluids. We study the long 

•^ ' term behavior of solutions and prove that the system possesses a pullback exponential 

attractor. In particular the regularity estimates we obtain depend on the initial data 
only through fixed powers of their norms and these powers are uniform with respect to 
the growth of the polynomial potential considered in the Cahn-Hilliard equation. 

< 

^ ■ 1. Introduction 

The modeling of multi-phase flows has been a fleld of intense mathematical research 
in the last years. In particular, its relevance for applications (see, e.g., [17, 18, 19, 21] 
and references therein) has motivated several different approaches, among which diffuse 
(^ . interface methods seem to combine both numerical efficiency and theoretical tractability. 
^ ! One of their instances is given by the so called model H, which was first proposed in [19] 
Q I and [22] and then rigorously derived in [17]. In the corresponding system, the two- 
^ ■ phase fiow is described by a (mean) velocity field u, which satisfies a Navier-Stokes type 
O ■ equation, and an order parameter field ijj, which represents the difference of the relative 
^^ ' concentrations of the two fluids and which solves a convective Cahn-Hilliard equation. 
More precisely, given a bounded and smooth domain fi C M^, assuming that the viscosity 
of the mixture is a constant z/ > 0, the two-dimensional model H reads as 



> 



X 



:l.l) 



dtu + u ■ Vu — z/Au = Vp -|- fiV^ + g, in D., 
V-u = 0, in fi, 

dfip + (u ■ V)ip = rriAfi, in Q, 
^/i = -eA^ + i/(^), infi. 



where /i is the so called chemical potential with constant mobility ?Ti ^ 0, g is a time- 
dependent bulk force, / is the derivative of a double-well potential F while p = ir + 

-\Vip\'^ + -F{^p) introduces the pressure vr as well as the positive parameter e rendering 

the interaction between the two phases. In particular, e is related with the small but not 
negligible thickness of the interface. 
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2 STEFANO BOSIA AND STEFANIA GATTI 

This system is usually complemented by homogeneous Dirichlet boundary conditions on 
the velocity field, no fiux boundary conditions on the order parameter field and chemical 
potential, namely, 

(1.2) u = 0, d^tp = 0, 9n/i = 0, on dQ. 

Being the problem non-autonomous, we specify the initial values at a given time r G M 
for the state variables, that is, 

u(r) = uo, ^(r) = ijQ, in Q. 

We recall that, in this model, the chemical potential of the binary mixture fi is given by 
the variational derivative of the free energy functional for the Cahn-Hilliard equation 



J'ii^) = J {^W\' + VFW) dx. 



where F{^p) is a suitable double- well potential characterizing the phase decomposition of 
the mixture. Since F' = /, the fourth equation in (1.1) follows. 

The Cahn-Hilliard model for spinodal decomposition and coarsening during quenching 
of alloys was first proposed in [8]. In this setting a thermodynamically consistent double- 
well potential F is naturally seen to be logarithmic (see [7] and references therein). How- 
ever, this singular form for the potential causes major difficulties in the numerical and 
theoretical study of the Cahn-Hilliard system so that in applications it is often replaced 
by a polynomial approximation like 

In this context, a possible approach to deal with the physically relevant case consists in 
suitably approximating the singular potential by polynomials of increasing order (see [11, 
12] for an application of this technique to a system closely related to the model H). Our 
paper deals with a polynomial potential F of arbitrary order p + 3, for p ^ 1 (in fact, the 
lower order case is much easier). 

From the mathematical viewpoint, system (1.1)-(1.2) has been firstly studied in [23] for 
Q = M^. Then, in the case of bounded domains, global existence results for both weak and 
strong solutions in the 2D case were obtained in [5] (see also [6]). More recently, the case 
of logarithmic potentials has been considered in [2] (see also [1]), where, in particular, 
the convergence of solutions to a single equilibrium has been established in absence of 
nongradient external forces. This issue has also been investigated in [24] for smooth 
potentials. A rather complete picture of the longtime behavior in the case A^ = 2 on a 
bounded domain can be found in [13]. In the case A^ = 3, existence of trajectory attractors 
has been demonstrated in [14] with time-dependent external forces. Many related models 
have also been extensively studied in recent years. Among others we recall non-Newtonian 
Ladyzhenskaya fiuids (see [16] for the 2D case and [3] for some partial results in the 3D 
case), non-local interactions (see [11, 12]) and chemically reacting mixtures (see [4] and 
references therein). 

Concerning the mathematical theory of infinite dimensional dynamical systems, pull- 
back exponential attractors represent a new instrument recently introduced in the litera- 
ture (see [20]). The related theory combines both the advantages of pullback attractors in 
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the non-autonomous case (see [9] for a comprehensive introduction to the theory of pull- 
back attractors) and of the exponential attractors in Banach spaces in their most general 
form known today (see [10]). 

The main result of this paper is the existence of a pullback exponential attractor for 
the system (1.1)- (1.2). As a byproduct, we derive several regularity estimates for the 
solutions to the Cahn-Hilliard-Navier-Stokes system. These have an interest of their own 
due to their uniform structure with respect to the growth of the double-well potential F. 
Indeed, if the potential /(VO is assumed to satisfy |/(^)| ^ Cd?/'!^^^ -|- 1), we are able to 
control the solutions only by suitable powers of the norms of the initial data independent 
of p. This is not an easy task since the computations repeatedly involve f{ip) and its 
derivatives, which are naturally estimated as 

1/(^)1. ^ Cm^.^' + 1) ^ C{\^pf^'^^\A^pf^'/^ + 1), 

carrying the polynomial character of F directly into play. This obstacle is circumvented 
by suitably handling the nonlinear terms so that the dependence on the "shape" of the 
potential is confined to the multiplicative constants appearing in our results. 

Therefore, these estimates can be seen as a preliminary step forwards an effective 
approximating procedure able to deal with the more physically relevant case given by the 
singular potential. 

The plan of the paper goes as follows. In Section 2 we introduce the functional setting 
required to study system (1.1) and the main results obtained in this work. After recalling 
the theory of pullback exponential attractors in Section 3, we first derive basic energy 
estimates (Section 4) and then higher order regularity estimates (Section 5): in particular, 
our results are uniform w.r.t. the shape of F in the sense made precise above. We then 
derive continuity results and time regularity for solutions in Sections 6 and 7. Finally, in 
Section 8 we are able to check all the assumptions of the abstract results of Section 3 in 
the case of system (1.1), concluding the proof of our main theorem. 

2. Functional setting and main results 

We will denote by fi a smooth bounded domain of M^. The spaces Lp{Q) will be 
the usual Banach spaces of p-integrable functions with p G [l,oo]. The Sobolev-Hilbert 
space, which consists of fc-differentiable functions in the sense of distributions with square 
integrable derivatives, will be denoted by H''{Q). We shall use the bold symbols L^(r2) 
and H*'(r2) for the corresponding spaces of vector valued functions. The space of functions 
belonging to H''{Q) and vanishing on the boundary will be denoted by Hq{Q). Norms in 
the Sobolev spaces H''{Q) will be denoted by || ■ ||^fc, whereas we will use the shorthand 
notation | ■ \p for the norm in L^{Q) spaces, 1 ^ p ^ oo. In order to study the velocity 
field u we introduce the usual framework of divergence-free distributions, i.e., 

V = {ct)eC^{Q;R^) I V-0 = O}. 

Then we consider its closure under suitable distributional norms 

We also introduce the Leray projector P: L^(r2) — > L^jy(r2) mapping every element of 
L^(r2) to its divergence-free part. Furthermore, we will indicate by HJi^(r2) the dual 
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space of HQ^j^(r2). In Hq{Q) and HQ^jy(r2) we will consider the following norms 
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Finally we will denote by (/, (?) both the scalar product in L'^{Vt) (or L^(r2)) and the 
duality pairing between H~^{VL) and -ffo(fi) (or their vector valued analogues), the exact 
meaning being clear from the context. 

Since the second equation in (1.1) together with the boundary condition imply that 
the bulk integral of the order parameter is preserved by the evolution, we need to suitably 
account for this feature. First of all we define the mean value of / over the domain Vt as 



'''-¥\L 



/dx. 



denoting by / the mean free part of /, that is, 

/ = /-(/)• 

Thus, up to a shift of the order parameter field, we can always assume that the mean 
of ip is zero at the initial time and, due to the conservation of mass enforced by the 
Neumann boundary conditions, this will remain true for all positive times. Then the 
order parameter will belong to subspaces of U'iyt) and H^iVt) consisting of functions with 
zero mean, defined as 

Ll.^iSl) = {v& L^iSl) I {v) = 0}, iJ(%(^) = {^ e H\Q) I (v) = 0}. 

Here we can use Poincare's inequality (and some of its variants) at several stages when 
estimating the Sobolev norms of ijj. Indeed, the boundary conditions and the above 
definitions imply that 



^ = 0, 9n^ = on dn, A^= d^i) = 0, d^Aip = on dVt. 

I n Jn Jdn 

Therefore, all the norms HV^Hiifj, J = 1, . . . , 4 are equivalent to the L^-norms of the deriva- 
tives of order j. Moreover, Korn's inequality holds. Thus we have 

Mm ~ |V^|2, Mm ~ |AV^|2, M\h^ ~ IVA7AI2, ||^||i^4 ~ |A2^|2. 

Finally, the functional spaces for the whole solution (u, tp) are 

Ho = LL(fi) X Hl,^in) n, = H^,(fi) X /7f,o)(fi), 

which arise naturally in the study of the process generated by the solution of system (1.1). 
We can now list the assumptions on the potential F{ip), starting with some hypotheses 
concerning its regularity and growth: 

(A.l) Fe C5(M). 

(A. 2) F{y) grows at most polynomially fast at infinity, namely 

\f"{y)\<Cf{l + \yn, 
for some positive constants p and Cf. 
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(A. 3) The potential is coercive, i.e. there exist positive real numbers q and c/ such that 

(2.1) F{y)^Cf{\yr'^-l) 
holds for all y G M. 

Recalling that the potential F appears in system (1.1) only through its derivative, without 
loss of generality we can further assume that 

(A. 4) the functional F is strictly positive, i.e. F{y) > 0, Vy G M. 
We now give some additional assumptions concerning the shape of the double-well 
potential F. 
(A. 5) F{y) is a quadratic perturbation of a regular convex function defined on the whole 
M, that is, 

F{y) = Fo{y)-ay'' + -fy + (3, 

where Fq G C^(]R) is convex and a G M is a positive constant. 
(A. 6) Up to a suitable choice of the constants /3 and 7 in Assumption (A. 5), the convex 
part of the potential Fq satisfies: 

Fo(0) = F^(0)=0. 

In order to obtain higher order estimates having uniform structure with respect to the 
growth of /, we will assume that F behaves as a polynomial at infinity. In particular we 
will suppose 

(A. 7) The relation q = p + 1 holds in Assumptions (A. 2) and (A. 3). Moreover, for any 
/c = 0, 1, 2, 3, 4 there exists a positive constant c^ such that 

(2.2) |/«(y)|^c4l + F(y)^^ + ^-^V(. + 3,]^ VyGM. 

Remark 2.1. In the following Assumption (A. 7) will be necessary in order to estimate the 
U norm of the derivatives of the potential F in terms of some L^ norm of the potential 
itself. In particular, when dealing with higher order estimates we will often use the 
immediate consequence of Assumption (A. 7) 

p+2-fc 

for any ip such that F{iIj) G L^{Q) and for some constant C depending on k and Q. 

Remark 2.2. Throughout the paper, we will always assume that p ^ 2, being much 
easier the case when the potential F{y) grows at most as y^ at infinity. In particular, 
all estimates in the following sections hold for p ^ 1, except for (6.6) below. However, 
note that a suitable estimate for this term can be produced also in the case p E [1,2) 
(cf. Remark 6.1) under the assumption 

/(»)(y)^C \/yeR. 

The case p = 1 is particularly relevant for applications since the polynomial potential 
F{y) = {y^ ~ 1)^) which is often used in numerical simulations, falls in this setting. 

Finally, we assume the non-autonomous forcing term (symbol) g appearing in equa- 
tion (1.1) satisfies the following conditions: 

(B.l)gGLL(M;LL(fi)). 
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(B.2) g G Ll,^^{~oo,t; L^^^ll^)), for any t G R, that is, 

Mg(t) = sup/ |g(s)|2ds < oo, \/teR. 
(B.3) There exist to ^ 11^ ^^^ q > 2 such that g G L^[q^(— oo, toi L^i^(r2)), namely, 

^g,g(^o) = sup / |g(s)|^ds < CX). 
r^to Jr-1 

In this paper we will prove the following main results: 

Theorem 2.1. Assume that g satisfies (B.l) and (B.3) and let Ug(t, r) : "Ho — ^ "^o ^e the 
solution operator for the system (1.1). Then there exists a family Aijj^ = {J^Ug^it): t ^ 
to} of nonempty compact subsets of T-ii, which is a pullback exponential attractor for 
system (1.1) (see Theorem 3.2 below) in the topology ofHi. 

Corollary 2.2. Under the same assumptions of Theorem 2.1, if, moreover, g satis- 
fies (B.3) uniformly for to G M (as in (B.2)j, then the process Ug{t,r): "Hi — )■ "Hi has a 

family Aiu^ = {-^{/g('^)- ^ ^ I^} of nonempty compact subsets ofHi, which is a pullback 
exponential attractor for system (1.1) (see Theorem 3.3 below) in the topology ofHi. 

3. Exponential pullback attractors 

In this section we briefly review the theory of exponential pullback attractors as de- 
veloped in [20]. Below, (if, | ■ |) and (V, || ■ ||) will be two Banach spaces such that V is 
compactly embedded in H. Both spaces have a metric structure therefore, given any two 
nonempty subsets Di, D2 of the metric space X = H,V, the Hausdorff semidistance and 
distance are well deflned respectively as 

distx(-Di,-D2) = sup inf ||fi-W2|U 

viGDi V2&D2 

and 

dist5^"^(Di, D2) = max{distx(/^i, ^^2), distx(/^2, ^i)}- 
As in the usual case of exponential attractors [10] the key point of the argument 
is the introduction of a set of mappings which enjoy a suitable "smoothing property". 
This is responsible for the exponential convergence of the trajectories of the system to 
an exponentially attracting flnite-dimensional compact set (an exponential attractor) as 
soon as the trajectories have entered a sufficiently small neighbourhood of the attractor 
itself. Let 5, -fT G M be positive constants and let i? be a bounded and closed subset of V . 
Deflne 'Bs^k{B) to be the class of mappings S: V —^V such that the smoothing property 
holds on a ^-neighbourhood (in V) of 5, i.e. 

S{Os{B)) c B 

and 

\\Svi — Sv2\\ ^ K\vi — V2\ for all fi,t'2 G Os{B), 

where Os{B) = {v E V \ inf^^g^ \\v — w\\ < 5} is a (5-neighbourhood of the set B in V. 

We introduce a suitable class of family of mappings, which are the abstract, discrete- 
time, dynamical system representation of the evolution equations we will be interested 
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in. In particular, let no G Z be fixed and consider the class Ud{V, uq) of all families 
U = {U{m, n) \ n,m ^ Z,n ^ m ^ uq} of mappings U{m, n): V -^ V such that 

(1) U{n,n) = Id for all n ^ uq; 

(2) U{m, k)U{k, n) = U{m, n) for any n ^ k ^ m ^ Uq. 

When dealing with puUback attractors, only the evolution of the system up to the "present" 
time Uq is of interest. The key question is how perturbations of the system in the past 
affect the present dynamic and what actually is the state of the system observed. This is 
the reason why elements belonging to the class Ud{V,nQ) are defined up to time uq and 
not necessarily beyond. 

We can say that element of the discrete time class UdiV, Tq) possess a discrete time 
puUback exponential attractor in the sense made precise by the following theorem 

Theorem 3.1 ([20, Theorem 2.1]). Let hq e Z, 6 > 0, K > and B C V be fixed with 
B hounded and closed in V . Then, there exist positive constants Ci, C2, e and a only 
depending on V, H, 6, K and B, such that, for each U G Wrf(V,no) satisfying 

U{n, n — 1) G Ss.xiB) for all n ^ no, 

there exists a family Aiu = {A^;7(n) | n ^ no}, of nonempty subsets ofV, which satisfies 

a) Aijj is positively invariant i.e. 

U{m, n)Aiu{n) C Aiu{m) for all n ^ m ^ Uq, 

b) A^(7(n) <Z B is a compact subset ofV, with finite fractal dimension estimated by 

log2 N^{Aiu{n), V) ^ Ci log2 — \- C2 for all < e < e and any n ^ no, 

where N^{M.u{n),V) is the minimal number of e- balls in V, which are necessary to 
cover Aiu{n), 

c) M.U attracts B exponentially in a puUback sense i.e. 

disty(f/(m, n)B, Mu{m)) ^ de-^^"-") for all n ^ m ^ uq, 

d) for every integer k ^ 

■Mu{n + k) = JUxi^uin) for all n ^ no, 
where TkU{m, n) = U{m + k,n + k). 

Remark 3.1. The results in [20] also include robustness of a discrete-time exponential 
puUback attractor w.r.t. to a suitable metric in the space of discrete time processes. For 
simplicity we do not mention all the pertinent details here. However, we recall that this 
result is important in deducing the analogue continuous-time theory and, in particular, 
in obtaining continuity in time of an exponential pullback attractor. 

Having in mind the more relevant continuous-time setting, we now introduce a suitable 
analogue of the class Ud{V, no). Let to G M be any time, and consider the class U{V, to) of 
all families U = {U{t, s) \ s,t ^^,3 ^t ^ to} of mappings U{t, s): V ^ V such that 

(1) U{s,s) = /dfor all s ^ to; 

(2) f/(t, r)f/(r, s) = U{t, s) for any s ^ r ^ t ^ to. 
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In this setting, a natural way to introduce a smoothing property is to consider only those 
families U G W(V,to) such that there exists a positive time span tq for which 

(3.1) Uit,t-To)eSs,KiB) 

holds for all t ^ to- Thanks to Theorem 3.1, for any t ^ to the family f/* G ^^(V, 0) given 

by 

f/*(m, n) = U(t + rriTQ, t + titq) for all n ^ ??7. ^ 

possesses a discrete time exponential puUback attractor. 

In order to obtain a satisfactory dynamical description of the system also in the 
continuous-time case, we will need some additional assumptions on the time regularity 
and continuous dependence of the family U G W(y, to). In particular we will assume 

(H.l) Continuity w.r.t. the forcing terms: there exist real positive constants Co, eo and 
7 such that eo ^ To and that for all t ^ to. To ^ r ^ 2ro, ^ s ^ eo and v G Os{B) 

||f/(t, t-r)v- U{t -s,t-r- s)v\\ ^ Co\s\^. 

(H.2) Past continuous dependence on initial data: there exists a positive constant Cb 
such that 

\\U{t,t - s)v - U{t,t - s)w\\ ^ Cb\\v-w\\ 

for all v,w & B and any t ^ to, ^ s ^ 2tq. 
(H.3) Time continuity of solutions: there exist positive constants Cq and 7' such that 
for all t ^ to. To ^ r ^ 2tq, ^ s ^ eo and f G -B 

||f/(t, t - r)v -U{t-s,t- r)v\\ ^ C'^lsp' . 
We can now state the main result on exponential pullback attractors 

Theorem 3.2 ([20, Theorem 2.2]). IfU G U{V,to) satisfies (3.1) and Assumption (H.2), 
with B G V bounded and closed in V , then the family Aiu = {Mu{t) I t ^ to}; defined 
by 

Mu{t) = U t/(t, t- s- To)Mut-s-To (0) for all t ^ to, 

sG[0,ro] 

satisfies 

a) U{t, t)Mu{t) C Mu{t) for allr ^t^ to, 

b) M.T-ru{t) = -M-uit — T) for allr '^ and any t ^ to, where T^T-U{t, s) = U{t — T, s — t), 

c) for all T ^ and any t ^ to 

disty(t/(t,t - T)B,Mu{t)) ^ Ce-"\ 

d) if, for any D (ZV bounded, there exists a time Sd ^ such that 

U{t,t- sd)D C B for all t^ to 
then 

(3.2) disty(f/(t, t - t)D, Muit)) ^ Ce^^^^e"^^ for all r ^ sd and any t ^ to. 
//, moreover, U also satisfies Assumptions (H.l) and (H.3), then 

e) M-uif) is a compact subset ofV, with finite fractal dimension, for all t ^ to, 
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f) for all ^ r ^ eo and any t ^ to; 

Aisi'^'^{Mu{t),Mu{t - r)) ^ C\r\\ 

Finally consider the case of processes on V. Let t/ be a family U = {U(t,s) | s,t G 
M, s ^ t} of mappings U{t, s): V -^ V such that 

(1) U{s,s) = Jdfor all seR; 

(2) f/(t, r)U{r, s) = U{t, s) for any s ^ r ^ t. 

This corresponds to a dynamical system defined not only up to the present time to, but also 
for positive times. Considering processes corresponds to investigating what the eventual 
fate of the system under scrutiny will be. Therefore, it is interesting to investigate the 
relation between this eventual fate and the present state of the system, which is itself, in 
a way of speaking, the outcome of an arbitrary long evolution. 

We will need the following additional assumption, which is a slight modification of (H.2) 

(H.4) Future continuous dependence on initial data: for any t > Iq and Di, D2 bounded 
subsets of V , there exists a positive constant L(t, Di,D2) such that 

\\U{t,to)v-U{t,to)w\\ ^ L{t, Di, D2)\\v - w\\ for all v E Di, w e D2. 

Theorem 3.3 ([20, Theorem 2.3]). Assume that U is a process on V and, for some 
to E M, the subfamily of U given by the operators U(t,s) when s ^ t ^ to satisfies (3.1) 
and Assumption (H.2), with B G V bounded and closed in V . Under these assumptions, 
the family Mu = {Mu{t) | t G M} defined by 

""^^ \uit,to)Mu{to) tft>to, 
where Aiu is the family given in Theorem 3.2, satisfies: 

a) U{t, T)Muir)^ Muit), for all r ^t, 

b) Air^^uit) = -Muit — t) for all r ^ and any t ^ to and 

^T-ruif) C Mu(t — t) for all T ^ and any t > to, 
where T^T-U{t, s) = U{t — t,s — r). 
If in addition (H.4) holds, then 

c) if, for any D G V bounded, there exists a positive time sd such that 

U{t, t — s)D G B for all s ^ sd and any t ^ to, 
then Aiu satisfies (3.2) for all t ^ to and 

disty(f/(t,t - T)D,Mu{t)) ^ Z(t,5,>i'c/(to))e°('^+*-*°^e-°" 
for all t > to and any T^S£, + t — tQ. 
Moreover, if U also satisfies Assumptions (H.l) and (H.3), then 

d) A^[/(t) is a compact subset ofV with finite fractal dimension for all t G M, 

e) for all ^ r ^ eQ and any t ^ to 

disf^"^ {Muit),Mu{t - r)) ^ C\r\\ 
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Remark 3.2. We recall that in [20] also explicit estimates on a fractal dimension of the 
pullback exponential attractor have been derived. For the sake of simplicity, we neglect 
them here. 



4. Existence results and basic energy estimate 

In this section we recall some basic energy estimates, which are obtained naturally 
when proving existence of solution to system (1.1). First of all, for the sake of simplicity, 
we set e = m = 1 and we write the definition of weak solution to system (1.1). 

Definition 4.1. Let zq = {uo,ipo) G L^i^(r2) x HLq^Q) and let r G M. Then a couple 
z = (u, ip) such that 

u G L'iT,T;Ul^,^mnH\T,T;ii2m 
tfj G L2(r, T; H^n)) n H\t, T; H~\n)) 

is called a weak solution to (1.1) if 

{^Mt),^r) + ((u(t) • V)u(t), v) + (z/Vu(t), Vv) = (^(t) V^(t), v) 
{dtm, 0) + (u(t) ■ Vm, 0) = - (V/i(t), V0) 
hold for a.e. t G [r, T], for all v G V and for all (p in C°°(f]), if 

fi{t) = /(^(t)) - A^(t) 
holds for a.e. t G [r, T] in H\n) with /i G ^^(r, T; i/^(fi)) and if 

limu(t)=uo inLdi^(fi), lim ?/^(t) = V^o in i^L)(^)- 

The well-posedness for problem (1.1)-(1.2) is justified in a suitable Galerkin scheme, 
thanks to the following a priori estimates and the subsequent Lemma 6.1(see e.g. [5, 13]). 

Theorem 4.1. Let assumptions (A.1)-(A.6) hold. Ifg satisfies (B.l) andz^ = (uo,'?/'o) ^ 
"Ho; then there exists a unique weak solution z(t) = (u(t), ?/^(t)) departing at time r from 
the initial datum zq. 

We now obtain the first basic energy estimates that will be the basis for the estimates 
of the following sections. 

Lemma 4.2. //g satisfies (B.l) and z{t) = {u{t),tlj{t)) is the solution departing at time 
T from the initial datum Zq = (uq, ipo) G "Ho? denoting by n{t) the corresponding chemical 
potential, there holds 

(4.1) \uit)\l + \Vm\l + 2|F(^(t))|i + ^ [z/||u(.)||^.^^^^(^) + |V/i(s)|2l ds 

^|uo|^ + |V7/;o|^ + 2|F(^o)|i + C / \g{s)\lds. 
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Besides, there hold 

(4.2) I* \V^{s)\lds + J' \F{^{s))\, ds + J' \Aij{s)\lds 

^C ("luol^ + |V^o|2 + 2\F{iJo)\i + j \g{s)\lds^ + C{t - t) 
as well as 

(4.3) I' \Ai;{s)\tds 

2 



^C ("luol^ + iVV^ols + 2|F(^o)|i + ^ |g(s)l2ds') 



+ C{t - t) ( |uo|^ + |VV^o|2 + 2|F(^o)|i + j Isismds 

Proof. In order to obtain our first (dissipative) a priori estimate, we multiply the first 
equation in (1.1) by u and the third by /i. Recalling the antisymmetric property of the 
convective term in the Navier Stokes equation and exploiting the useful vector identity 

= ^^(|V^|^ + 2|F(V;)|i), 

we obtain 

(4.4) ~ (|u|^ + IV^I^ + 2\Fm,) + HVu|^ + |V/i|^ = (g, u) . 

Recalling Poincare inequality for u, integrating this formula with respect to time, we 
deduce (4.1). 

We now have to "complete the norms" on the left hand side of (4.4). From the definition 
of the chemical potential /i (i.e. from the fourth equation in (1.1)) we have 

Since, by assumption, ip is mean free, we also deduce 

{fi,^) = {fi-{fi),^)^^\Wfi\l + C\^\l, 

where C is a constant, which only depends on the domain Q. From assumption (A. 5) on 
the potential F we further deduce 

(/(^),^) = (/o(^),V^)-2a|V;|^, 

being /q = Fq. Taking into account the convexity of Fq we can also bound the right hand 
side of this identity from below: 

{foiij),i;)^\Foii;)-Fomi- 
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Putting the last four estimates together and recalhng Assumption (A. 6), we obtain 

= |V^|^ + |F(^)|i-a|^|2. 
Therefore we get 

\Vij\l + \Fmi ^ |V/i|^ + C\^\l^ \V^^\l + 5\i,\lXl + C, 

where g is a positive real number, 5 is a (small) positive constant, which will be determined 
later, and C is a positive constant, which depends only on the domain Vt and is independent 
of the exponent q as soon as g ^ g > 0. 

Remark 4.1. We observe that under assumption (A. 7) we immediately have g ^ 2 so that 
in our case the constant C really depends only on Vt. 

By adding this last estimate and (4.4) together, choosing 5 small enough, we finally 
deduce the basic energy estimate for system (1. !)-(!. 2) 

(4.5) ^ (|u|^ + IV^I^ + 2\Fmi) +C{\Vm\1 + \V^\l + 2|F(^)|i + | V/.|^) 

^^(l + |g|^)- 

Integrating with respect to time from r to t, we then obtain the first part of estimate (4.2). 
Noticing that 

(V/x,V^) = -(/i,A^) = |A^|2-(/(^),Az^) = |A^|2+(/'(V;)VV^,V^) ^ |AV;|^-2a|V^|L 

we have 



(4.6) |AV'|^^|V/i|2|V^|2 + 2a|V^|2 



2' 



which, integrated in time, on account of the above estimate (4.1), gives the second part 
of estimate (4.2). 

In order to prove (4.3), we square (4.6), obtaining 

\Ai,\l ^ C{\V^^\l\Vi^\l + IV^I^) ^ C{\V^^\l\Vi,\l + IV^I^IA^I^). 

By an integration in time, in view of (4.1) and (4.2) we accomplish our purpose. D 

Corollary 4.3. // g satisfies (B.l) and (B.2) and z{t) = {u{t),ilj(t)) is the solution 
departing at time r from the initial datum Zq = (uq^iPq) G Hq, the following dissipative 
estimate holds 

(4.7) |u(t)|^ + |V^(t)|^ + 2|F(^(t))|i 

^ (|uo|^ + I V^ol2 + 2|F(V;o)|i) e"^(*--) +0(1 + M^{t)) , Vt ^ r. 
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Proof. The dissipative estimate easily follows from the basic energy estimate (4.5) using 
Poincare's and Gronwall's inequalities as well as the known estimate 
rt 



e 


-"t 


.^1g(.)|^d. 


G 


oo 
n=0" 


/t-(n+l) 


c 


'sup / 


\s{s)\lds, 
-1 



g(s)|2ds^e-^*ye^(*-")sup T \g{s)\lds 



which holds for g G Lf,iQ^(-oo, t; L2(n)). D 

Remark 4.2. A bound on VA?/' in L^(r, T; L^(r2)) can also be easily deduced by computing 
the L^ norm of the gradient of the equation for the chemical potential fi in (1.1) thus 
leading to the regularity of the order parameter field required by Definition 4.1. However, 
this estimate cannot be easily made uniform with respect to the shape of the potential F 
(and in particular with respect to the growth exponent p). 

Remark 4.3. From the above computations we deduce the following regularity for weak 
solutions of system (1.1) 

u G L^iT,T; LL(fi)) H L^ir^T; lil^,^m 
^ G L°°(r, T; H\Q)) n L\t, T; H\Q)) 

V/iGL2(r,r;L2(fi)) 
for any T G M, T > r. 

5. Higher regularity estimates 

In order to obtain estimates having uniform structure with respect to the growth ex- 
ponent of /, we henceforth assume that F satisfies (A. 7). Although all exponents and 
norms that appear in this and in the following sections are independent of p, the general 
constant C will quickly become larger as p grows. 

In particular, Assumption (A. 7) and Lemma 4.2 imply 



(5.1) \f^'\m)\^^ ^ci\Fim)\i + i) 

^ C ("luol^ + \VMl + 2|F(V;o)|i + ^ Us)\lds + l\ , 

being (u(t),'?/^(t)) the solution to (1.1)-(1.2) departing from (uq, ^o) ^ "^o ^^ time r. 

The goal of this section is to improve "by one order" the basic regularity result already 
obtained. In particular, under suitable assumptions, we will get to u G L°°(Hj jjy(r2)) fl 
L2(H2d;^(fi)) and %Ij G L'^{H'^{VL)) n L^{H^{n)). This will be achieved in several steps 
gaining before spatial regularity for f{ijj) and /i and later time regularity as well: first in 
Lemma 5.1 we will deduce f{ip) G L'^{L^{n)) and fi G L'^{L'^{n)); then f^tp) G L'^{L'i{n)) 
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and A^ G LP'{L'^{VL)) for any g ^ 1, as shown in Lemma 5.2; this will give /i G L°°{L'^{VL)) 
(cf. Lemma 5.3) and the final result (see Lemma 5.4). 

Notation. In order to simplify notation, we will denote by At^r the quantity 

At,r = 1 + |uo|^ + |V^o|2 + 2|F(7/;o)|i + / |g(s)|^ ds, 

which depends only on the initial data Uq, ipo, on the forcing term g and on the times t 
and r. Besides, C stands for a generic positive constant depending only on Q and possibly 
on p and is allowed to vary even in the same line. 

Lemma 5.1. If z(t) = (u(t),'?/;(t)) is the solution departing at time r from the initial 
datum Zq = {uq,iIjq) G Hq, denoting by /x(i) the corresponding chemical potential, there 
holds 

(5.2) j\\f{^{s)\l + \^^{s)\l) As ^ CAl + C{t - r)A,. 

for any t '^ t, t eM.. 

Proof. Remark 2.1 allows to bound the mean value of fiip) as 

\{f{m^c\fm^^c{i + \Fm,), 

p+2 

for some C > depending on p only through the constant cq in Assumption (A. 7). 
Recalling the equation defining the chemical potential in (1.1) and estimate (4.6), we 
further deduce 

\fW - {fm I2 ^ 2|^ - (/i) \l + 2| A^|2 ^ C| V^l^ + 4a| V^l^. 

Therefore, we gain full control on the L^-norm of f{ip), bounding its time integral as 

\fWs))\lds 
^C f \Vfiis)\lds + Aa f \\/i>is)\lds + C f il + \F{^{s))\,yds 

where we used (4.1) and (4.2) from Lemma 4.2. The second part of estimate (5.2) follows 
from 

|/i|^^2(|AV;|^ + 1/(^^)1^), 
the above bound and (4.2). D 

As announced before, the integrability of f{ip) can be further improved in two steps. 

Lemma 5.2. If z(t) = {u{t),tp{t)) is the solution departing at time r from the initial 
datum Zq = (uq, ^0) £ ^0, then, for any b > there exists Cb > such that 

{\fWs))\h2 + \^^{s)\U2) ds ^ CbAl + a(t - r)A,., 
meaning that fiip), Atp G L'^{T,t] L^~^'^{Q)), for any b > 0. 
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Remark 5.1. We note that this estimate extends to singular functional /, without appeal- 
ing to approximation arguments as in [2] but with the same order of controL 

Proof. Multiplying the equation for the chemical potential by f{ip)\f{ip)\^ and integrating 
over Q, we have 

(5.3) (/i,/wi/wr> = i/wiK-(A'0,/wi/wi'>- 

We now exploit assumption (A. 5) on F, proving after an integration by parts 

- (A^;, fmfm') = ib+i) (/'w, i/wnv^p) ^ -2aib + 1) {\fm\ iv^p> . 

Replacement of this estimate in (5.3) above leads to 

1/(^)1^:^ ^ 2a{b + 1) {\fm\ ivv^o + (i/ii, \fwD . 

Holder's and Young's inequalities then provide 
Recalling that from standard interpolation the inequality 

(5.4) l^l^:^^a|/^l^ll/^ll^i, 

holds, we end up with 

(5.5) |/(^)lK^^a(|V^|^|A^|^ + H^||/i||^.)- 
A further application of Young's inequality then gives 

Finally, integration with respect to time, leads to 

\fms))\l2ds ^ C,J^ {\V^is)\l + \A^is)\l + IM^)!^ + |V/i(.)|^) d., 

thus Lemmata 4.2 and 5.1 provide the first part of the desired estimate. To complete our 
argument, it is sufficient to exploit the equation for the chemical potential /i and this last 
estimate together with Lemmata 4.2 and 5.1: 

|A^|,%2 ^ 2 (|;.|^^2 + 1/(^^)1^+2) <H\f^\l + |V/i|^ + 1/(^)1^+2) • n 

Remark 5.2. Provided that F satisfies (A. 7), estimates (4.1) and (2.2) entail f{ip) G 
L°°(r,t;L(P+3)/(P+2)((])) (cf. Remark 2.1). Besides, the above Lemma 5.2 implies f{^) G 
L^{t, t] L^^^{Q)) for any b > 0. Being {p + 3)/(j> + 2) > 1, by the interpolation inequality 

(5.6) \hU^myh\l-l where «=^i^;j^M_, 
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we deduce that, when 4(1 — 6) = 4—- ; = 2, that is, b = 1 + p, 

2{bp + 2b + p + 1) 

t 



^CAl + C(t - T)Al. 

On account of (2.2), Lemmata 4.2 and 5.2, it thus follows f{tp) G L^(r, t; L^(r2)) and, in 
particular, 

(5.7) \\fm\lHr,t;LHn)) ^ CAl + C(t - r)^,. 

Lemma 5.3. If z(t) = {u{t),ip{t)) is the solution departing at time r from any initial 
datum zq = (uo,'(/^o) ^ ^i so that /ig = fi'ipo) ~ ^4'o ^ -^^(^); ^^c'^ there exists C > 
depending only on p such that the chemical potential fi is bounded in L°°{t,T] L'^{Q)) fl 
L'^{t,T; H^{n)) for all T > r and there hold 

Ht)\l ^ C (IMr)l^ + Al + (t - r)A,) e^(^tH-(*-)4.) 
and 

1^ |A/.(.)|^ ds^C {Hr)\l + Al + {t- r)A,0 (4. + (^ - ^)4.) e^(^-+(*--)^-). 

Proof. The evolution of the chemical potential yU is governed by 

(5.8) dtfi = f'WAfi - f{ij){u ■ V)^ - AV + A((u ■ V)^), 

as can be seen by formally differentiating with respect to time the last equation in (1.1) 
and by taking into account the third one. The product of this equality by fi gives rise to 
three terms from the right hand side: in order to exploit the lower bound on /', the first 
one can be written as 

{f'WAf,,^,) 

= - (/'(v;)v/i, v/i) - (r(^)v^,/iV/i) 
= - (/'(v^) v/i, v/i) - ^ (r (^)vv^, v(/i2)> 

= - (/' (^) v/i, v/i) + 1 (r (^) AV', /i^> + 1 {rm i vvi' , ^^'') . 

Thanks to the incompressibility condition, the second one reads as 

-(/'(^)(u-V)^,/i) = (/(^),u-V/i). 

Finally, noticing that the third equation in (1.1) and the boundary conditions (1.2) imply 
9n(A/i — u ■ Vip) = dndtip = on dQ, the last term is 

- (AV - A(u ■ V^), /i> = -|A/i|2 + (u ■ V^P, A/i) . 
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These computations lead us to 

(5.9) ^^|/i|^ + |A/i|2 

= - (/'(^) v/z, v/i) + 1 (r (v^) A^, /i^> + 1 (r (v^) I v^i^ , /z^> + (/(^)u, v/i) 

+ ((u-V)^,A/i). 

By Assumption (A. 5) on the potential F, the first term on the right hand side of iden- 
tity (5.9) is easily controlled by 2a|V/i|2 while the last one can be bounded by 

|((u-V)V^,A/i)| 

^|u|4|VV'|4|A/i|2 

^c|ui;/^ivui;/^ivv'i;/=^iA^i;/^iA/i|2 

^^|A/i|2 + C|u|2|Vu|2|V^|2|A^|2. 
o 

We are left to consider the other terms in (5.9). Having in mind (5.1), we prove 

\\{f"WA^,fi')\ 



■2' 



2 

m 



^C|r(^)U|AV;U3(|/i|^ + |/i|2|A/i|2) 
p 

^l\Afi\i + c [\fmi^ + 1) (iA^i;+3 + i)i;.|^. 



Analogously we have 



There also holds 



i|(n^)ivv^p,/i^>l 
p-i 

^C\rm^\A^|J\li\^^\l + \^^\2\A^^\2) 

p-l 

<l|A^i^ + ciirmu + mA^pit + i)i/i|^. 

O p-1 



|(/(^),u-V^)| 

^|/(^)|p±3|u|2(p+3)|V/i|2(p+3) 

p+2 

^C|/(^)U|Vu|2||/x||h2 

P+2 

^l\Af,\l + C\Vf,\l + C\fmU\^^\l 

O P+2 
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Collecting the above estimates and recalling (5.1), we have 

(5.10) ^^\^\l + \A^,\l<:h\fi\l + g, 

where 

/. = C(1 + |F(^)|0^(1 + |A^|^ + |A^|J^3) 
g = C\Vfi\l + CluyVuyVVblAV'b + C(l + \F{i,)\if\Vu\l 
In view of (4.3), Lemmata 4.2 and 5.2, h and g are integrable quantities. Indeed, we have 



^ h{s) ds^C {Al + {t- T)Al) , 



By Gronwall's lemma we further deduce 



Ht)\l ^ {Hr)\l + I g{s)ds) exp ( I h{s)ds 



that is, 

W)\l ^ C {\^^{r)\l + Al + (t - r)A,) e^(^t^-(*-)4.). 
Moreover, integrating (5.10), from the estimates above we also deduce 



^Hr)\l + l\{sMs)\lds + l\{s)ds 

^C {Hr)\l + At + {t- r)A,0 {Al + {t - r)Al) e^(<^+(*-)^t). Q 

Remark 5.3. The above Lemma has several consequences. First of all, from the third 
equation of (1.1) we easily obtain l^t^H ^ C(|A/i|2 + |u|2|Vu|2|VV'|2|A'?/'|2). Thus, Lem- 
mata 4.2 and 5.3 yield 

(5.11) j \dMs)\lds 

^ C {Ht)\1 + Al + (t - r)A,0 {Al + (t - T)Al) e^K.(-^K.). 
Besides, by (5.5) with 6 = 2, we have 

\fmt^C{\V,p\l\Aij\l+\f^\t+\fi\l\Vii\l), 
hence Lemmata 5.1 and 5.3 entail f{tp) G L^(r, t; L^(r2)), with 
ft 



^C{Ai^ + {t-r)Al^) _ 



(5.12) / |/(^(s))|^d. 

^ C {Ht)\1 + Al + {t- r)A,0 {Al + {t- r)A,0 
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Remark 5.4. Actually, even more uniform estimates can be deduced from the above Lem- 
mata. For example, from (5.5), using Ladyzhenskaja inequality and interpolation esti- 
mates, we deduce 

\f ml ^c{\\/^ij\t\A^p\t + \f^\i + \fi\t\\/ fill) 

^C{\V^|J\t\A^|J\t + \^^\l + \^^\t\A^^\l), 
i.e. /(V^) G L^{T,t;L^{Q)). In particular, 

(5.13) r \fiHs))\lds ^C f (|V^(.)|^|AV;(.)|^ + |M^)I^ + IM^)l^|AM^)l2) ds 
Jt-i Jt-i 



^C{W-l)\l + Al_,YA, 



t,t-i 



e 



Thanks to our assumptions on / and to the previous results, we can now obtain esti- 
mates on the higher norms of the solution, which have uniform structure w.r.t. the shape 
of the potential. Dependence on the growth of the potential F is limited to the constants 
C, which appear in the estimate. 

Lemma 5.4. Given any initial datum zq = (uo,'?/'o) ^ "^i so that Ho = f{ipo) — Aipo G 
L^(il), the solution departing at time r from Zq satisfies 

(5.14) ||z(t)||^^ 

^C{\\z4l, + If^oll + Al + (t - r) AO {Al + (t - r)Al) e^(^t.+(-^Mt) 

for some constant C > depending on the exponent p and on the domain Q but indepen- 
dent on the initial data. Moreover, 

{\A',Pis)\l + \Anis)\l)ds 

< C (llzoll?,, + l/iol^ + Al + {t- r)A,0 {Al + {t- r)Al) e<^t.+(*-^)4.). 
Proof. In this proof, we will exploit the following inequality 

\f'iy)\ + \ny)\^Cil + \f{y)\), VyGM, 

which can be easily obtained from Assumptions (A. 3), (A. 4), (A. 5), (A. 6) and (A. 7) by 
means of the Young's inequality. In particular, we will take advantage of its straightfor- 
ward consequence 

if'wu+irm^^cii + ifm,). 

Adding together the product of the first equation in (1.1) by 2Au = — 2PAu and of the 
third one by 2A'^ip, we obtain 

(5.15) 

^^{\Vu\l + \Ai,\l)+2u\Au\l + 2\A'^P\l 
=2 (g, Au) - 2 ((u ■ V)u, Au) + 2 (/iV^, Au) - 2 ((u ■ V)i/j, A^^) + 2 (A/(^), A^i/j) . 
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The first two terms arising from Navier-Stokes equations can be dealt witli by writing 

|-2((u- V)u,Au) + 2(g,Au)| 

^C(|u|;/^|Vu|2|Au|f + |g|2|Au|2) 

^^\Au\l + C\u\'\Vu\t + C\g\l. 

Since {fiVip, Au) = — {Atp'Vtp, Au), from the Agmon's inequahty we easily have 

|2(/iV^,Au)| 
^2|AV'|oo|VV'|2|Au|2 

^ciAv;i;/^iAvi;/'iv^i2iAui2 

Then, by Ladyzhenskaja inequality and standard estimates 

^2|U|4|V^|4|A2^|2 

^c|ui5/^ivui;/^ivv^i;/^iA^i;/^iAVi2 

^^\A'i,\l + C\u\l\Vn\l + C\V^l^\l\A^P\l. 
We are left to consider the last term in (5.15), for which we exploit 

namely, 

2\{Af{^p),A'^p)\ 
= |2 (/'(V^)A^, A^^) + 2 (r(^)|V^r, A^^) I 
<:C{1 + |/(V;)|4)(|A^|4 + |V^U|V^|4)|A2^|2 

^C(l + |/(^)|4)(|A^|;/^|VAV;|;/^ + \Vij\T\VAi:\f\Vi^&'\A^\l/')\A'ij\, 
^C(l + |/(^)|4)(|V^|f lA^^i;/^ + \V^\T\A'ijtt)\A'^\, 

We finally deduce the differential inequality 

d 
(5.16) 



^(Vul^. 


\Ai,\l) 


+ u Aul+ AVI2 


^Cu\lVu\l 






+ C{\ul 


Vu|2 + 


(1+ fwm^^i 



(|V^|^ + l)|A^|^ + |g|^). 



l + V ' ^'^"''l^ ' ^^'^-''l^ ' I"l2^ 



Introducing 

h = C\u\l\Vu\l 

g = C (|u|2|Vu|2 + (1 + |/(z^)|l)(|V^|^ + 1) + (IV^I^ + l)|Az^|^ + |g|^) 
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the above differential inequality reads as 

(5.17) ^^Mk^hMl^+g. 

Thus Gronwall's lemma gives 

l|z(t)||?,, ^ (||z(r)||^^ + ^^(s)ds)e/>W'^^ 

where, by (5.12) and Lemma 4.2 we have 



h{s) ds ^ CA^ 



t,T 



j\{s) ds^C {\^,{T)\l + Al + (t - r)A,0 {Al + (t - r)0 e^(4.+(*-K.) 
and the estimate (5.14). Moreover, integrating (5.16), we also have 

£(|A^z^(.)|^ + |Au(s)|^)d. 

^ C (llzoll?,, + |^(r)|^ + Al + {t- r)A,0 (4. + (^ " -)4.) e^(^-^(*--)^-). □ 

In the case of regular initial data, i.e. Zq G "Hi, we thus have the sought higher regularity 
for solutions: 

u e L-(r,T;H^,(n)) nL2(r,T;H^,,i,((])) 
^P e L'^ir, T; H\fl)) n L\t, T; H\n)) 
dt^ljeL''{T,T;L\n)) 

Corollary 5.5. Given any symbol g satisfying (B.l) and (B.2) and any to G M, t/iere 
exists a positive constant Cg(to) snc/i t/iat, /or any bounded set D C Tio, there exists 
T = T{\D\) > depending only on \D\ = max{l, sup^g/j ||z||-^„} such that 

Mt)\\n, + |/i(t)|2 ^ Cg(to), t ^ to, r ^ t - 4 - T(|D|), 

where z(t) zs t/ie solution to the problem with symbol g, departing at time r from the initial 
datum Zq G D, and fi{t) is the corresponding chemical potential. Besides, the following 
integral estimates hold true 

f |/(^(a)|^ds^g(Mg(to)) 
Jt-i 

fort^to,Ti^t-A-T{\D\), and, for t ^ to, r ^ t - 5 - T{\D\), 

f (lAu(.)l^ + |AV(5)|^) ds ^ Q(Mg(to)), 

for some nonnegative increasing function Q depending on p only through a multiplicative 
constant. 
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Proof. In order to prove the claim, we divide our argument in several steps. At each step, 
thanks to the estimates of the previous sections, we will be able to improve the regularity 
of the solution of (1.1)-(1.2) by letting the system evolve for a time interval sufficiently 
large (with the only exception of the initial step, however, all time steps will be taken 
equal to 1). 

Firstly, thanks to Corollary 4.3, for any to G K, any symbol g as above and any bounded 
set D C Tio, there exists T = T{\D\) > such that 

\\z{t)\\l^^ + 2\F{m)\i ^ 1 + C(l + Mg(to)) ^ C(l + Mg(to)), t ^ to, r^t-T, 
for any zq G D. Moreover integrating (4.5) and (4.6), for t ^ to 

/ {Hs)\\l^ + \V^^is)\l) ds ^ C(l + Mg(to)) + ||z(t-l)||?,.^ + 2|F(VKt-l))|i, r ^ t-1, 
so that, provided that r ^ t — 1 — T, we deduce 

(5.18) / (||z(s)||^^ + |V/i(s)l2)ds^C(l + Mg(to)) 

as well as 

(5.19) A,, ^ C(l + Mg(to)), r + T^t-l^s^t^to. 
It thus follows from Lemmata 4.2, 5.1 and 5.2 

r (lA^(.)l^ + \A^is)\l_,, + Hs)\l) ds ^ CAl_, ^ C(l + Mg2(to)), r ^ t - 1 - T, 
Jt-i 

allowing to prove that the functions h and g in the differential inequality (5.10), for 
r + T<t — l<t^to, satisfy 



^C(l + M,^(to)). 



r h{s)ds ^ CAl_, ^ C{1 + Mg^(to)), r 9{s)ds ^ CAl_, 

Jt~i Jt-i 

Hence, by the Uniform Gronwall's lemma and (5.19), it follows 

(5.20) |/i(t)|2^C(l + M|(to))e^(^+*'-*M, r + T + 2^t^to 

as well as 

/ \Afxis)\lds ^ C (1 + Mg^(to)) e^(^+^'-'(*°)), r + T + 3 ^ t ^ to. 
Jt-i 

Then the first claimed integral estimate follows from (5.13) and (5.20), while the functions 
h and g in (5.17) satisfy 

/ h{s)ds^C{l + Ml{to)) I (7(s)ds^C(l + M|(to))e^(i+*^-'M, 

provided that r + T + 3^t^to. Therefore, applying the Uniform Gronwall's lemma 
to (5.17) we deduce 



,^C2(to), r + T + 4^t^to, 

where 

C^(to)=C(l + M|(to))e^(i+*^^'M. 
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Finally, provided that r ^ t — 5 — T, integrating (5.16) over {t — l,t) we obtain 

r (|Au(s)|2 + \A'ij{s)\l) ds<:C{l + Mf (to)) e^(^+*^-'(*»)). D 

Jt-i 

6. Continuous defendence 

In this section we obtain continuous dependence estimates of the solutions w.r.t. initial 
data and forcing terms (see Lemmata 6.1 and 6.2). 

In order to address the first issue and to simplify notation, throughout the section 
we indicate by symbols with no subscripts the difference between quantities denoted by 
subscripts 1 and 2, i.e., 

f = fl- /2. 

From (1.1)-(1.2) we easily see that the difference between two solutions zi = (ui, ipi) and 
Z2 = (u2, ^2) satisfies the system 

' dtU + (u ■ V)ui + (U2 ■ V)u - Z/Au = fXiWl/jl - IX2^'4^2 + g 

V-u = 

(6.1) < atV' + u- V^i + U2 ■ VV' = A/i in Vt, 

fx = -A^P + /(V^i) - /(V^2) 
^^(r) =^01 -^02, u(r) = Uoi-Uo2 
with the boundary conditions 

u = 0, dni/j = 0, dnfi = 0, on dQ. 

Lemma 6.1. Let Zoi,Zq2 G "Hq be any pair of initial data and let gi,g2 ^ L^^ci^'^^^^'i^'divi^)) 
be any pair of symbols. Then there exists a constant C such that, if Zi{t), i = 1,2, are the 
solutions of (1.1)-(1.2) with initial data Zqj at time r and symbol gj, then the following 
estimates hold 



2 ^ / ll„ l|2 , / |„^„m2 j„ \ ^C" Mt,r +(i-T)ylt, 



-6 . ,, S-; 5\ 



and 



Hmiio^ [\\Mno + J Hs)\idsje 

l\\\/uis)\l + \VA^{s)\l)ds 
^C (^||zo||?,„ + I' \gis)\lds^ (a:^ + (t - r)^;') e 



where At^r is defined by At^^ = ^t,T(i) + ^t,T(2) ^'^^ ^t,T(i) ^^ the quantity corresponding to 
the initial datum Zj and the forcing term gj, for i = 1, 2. 

Proof. Recalling that, thanks to the incompressibility assumption on the velocity fields Uj, 
for i = 1, 2, we have 

= - {A^PiV^p^, u,) + {f(^p,)V^p^, u^) = - (A^.V^^, u^) + (VF(V'»), u,) 
= - {AtpiVipi, Ui) . 
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Moreover, the product of the first equation in (6.1) by 2u gives 

(6.2) — |u|^ + 2i/|Vu|^ + 2 ((u ■ V)ui, u) = -2 (A^Vz/^i, u) - 2 (A^sV^, u) + 2 (g, u) . 

Noticing that 

(A/i, A-ip) 
= {A[-A^p + f{^p^)-f{^P2)],A^P) 
= 1 VA^I^ + {f'i^P,)A^p, A^) + {[f^ - /'(^s)] AV^2, A^) 

+ {rMVij{Vij, + V^2),Aij) + {[f'M - r(^2)]|V^2p, AV> 
^|VA^|2 - 2a\A^P\l + ([/' (^0 - f'M]^iJ2, M) 

+ (r(^i)VV^(Vz^i + VV^2), A^) + ([/"(^i) - r(7A2)]|V^2|', At/;) , 

the product of the third equation in (6.1) by —2Ail) yields 

^1 VV^I^ + 2| VA^|2 - 2 ((u ■ V)V^i, A^;) - 2 ((u2 ■ V)'iP, A^) 
^4a\A^\l - 2 ([/'(^O - /'(^2)]A^2, A^) 

- 2 (/" (V;i)V^(V^i + ^2), A^;) - 2 ([/" (V^i) - r(^2)] I V^2r, A^) . 

Adding this last inequality to (6.2) we have 



d 



(6.3) - {\u\l + \Vij\l) + 2i.|Vu|^ + 2|VA^| 



dt 



^Aa\A^\l - 2 ([/'(V^i) - /'(V^2)]AV;2, A^) - 2 (/"(V;i)VV^(V7/;i + V^2), A^) 

- 2 ([/" (^i) - r (^2)] I VV^2|', A^;) + 2 ((U2 ■ V)^, A^) 

- 2 ((u ■ V)^, AV^2) - 2 ((u ■ V)ui, u) + 2 (g, u) . 

Since \A%Ij\2 ^ |V'?/'|2|VA'?/;|2, we immediately obtain 

4a|A^|2 ^ 4a|V^|2|VAV'|2 ^ ^IVA^I^ + CIWI'- 

5 

Before dealing with the terms arising from the double-well potential F, we introduce 
some useful estimates for terms of the form \f^^\il>i) — f^''' {"4^2)1 ■ Indeed from the funda- 
mental theorem of calculus and by assumptions (A. 7) and (A. 5) we have 






r>l/)2 /•1p2 



J1JJ2 J^i J ^1 

^C(l+F(^l)''' + '"'V(P + 3)+^(^2)(P+l-'=)/(p + 3))|^^_^^|_ 
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Moreover, we can estimate suitable L'^ norms of this difference as follows 



P+3 



(6.4) \f'Ki^i)-rK4^2)\^^ 

p-\-2 — k 
/ r X (P + 1 - fc)/(p + 3) / n N 1/(P + 3) 

^Cij (l + F(7Ai) + F(7A2))dxj ij |z^i-7/.2r'dxj 

^c (i + \F{^,)tr'~''"''''' + \F{i,,)tr-''"'^^''') i^u- 

In light of this bound, we can now resume the estimation of the terms in the right hand 
side of (6.3). We have 

|2((/'(V'i)-f(^2))A^2,A^)| 

^2|/'(^l) - /'(^2)|£±3|A^2|p+3|A^|,+3 
p+1 

^C (l + \F{^,)\t^''' + \F{^,)\t^'') iVV^blAv^^lp+slVA^b 
^ JlVAV^I^ + C (l + |F(^Olf ^"'' + \Fm\t^'''y |A^2|5-,3|VVli 
Analogously we also obtain 

l-2((r(^i)-r(V'2))ivv'2r,A^>l 

^2ir(v^i) - r(v^2)uiv^2ij+3iA^ip+3 

p 

<c(i+iF(^oir""/'^'^^^+ii^(^2)ir-"/'^"^')iv^yA^2i^ivA^i, 
^iivA^i^ + c (i + iF(v^oir '^'^"^' + ii^(^2)ir""/'^"^')' lA^^i^ivv^ii 



The last term involving the potential F and its derivatives can be dealt with in a similar 
way. Again the assumptions on F provide 

|-2(r(^i)VV'(VV'i + VV^2),AV')| 
^2|/"(^i)|p+3|V^|p+3 (iV^AiUa + |V7/;2|p+3) \M\v+z 

V 

^C (l + \Fi^,)\f'"-''^ |A^|2 (|AV;i|2 + |A^2|2) |VAV;|2 

^C (l + |F(V;i)|f ''"^^) IWIJ/^VAV^IJ/^ (lA^il^ + |AV;2|2) iVAz^l^ 

^^IVA^I^ + C (1 + |F(V;i)|t) (lA^il^ + \Aij,\t) \V^\l 

We also bound the last four terms on the right hand side of (6.3), which arise from the lin- 
ear momentum equation. Thanks to Agmon's inequality, we have | VV^I^ ^ C"! VV^|2| VA'?/^|2, 



26 STEFANO BOSIA AND STEFANIA GATTI 

and thus 

|2((u2-V)V^,AV;)-2((u-V)^,A^2)| 

^2|U2|2|V^|oo|AV'|2 + 2|u|2|VV'|oo|AV^2|2 
^C|U2|2|V^|2|VAV'|2 + C|U|2|VV^|;/'|VA^|;/'|AV;2|2 

^^\VAi,\l + ^|Vu|2 + C(|u2|^ + |AV'2|^)|V^|^ 

holds. Finally, by Ladyzhenskaja inequality, it follows 

|-2((u-V)ui,u) + 2(g,u)| 

^C|u|2|VUi|2 + C|g|2|Vu|2 
^C|u|2|Vu|2|Vui|2 + C|g|2|Vu|2 

^^|Vu|2 + C|Vui|2|u|^ + C|g|i 

Replacing the above estimates in (6.3), we see that ||z(t)|||^ = |u(t)|2 + |V^(t)|2 satisfies 

d 

dt' 
where h is given by 

/i=C(l + |Vui|2 + |Vu2|^) 

+ C (1 + |F(^i)li + inMt) (lAz^il^ + |A^2|^ + |AV^2|J+3) • 

By the results of the previous section we deduce 



(6.5) -||z||^^ + z/|Vu|^ + |VA^|2 ^ C\g\l + h\\z\\l^, 



f h{s) ds^C (a,/ + {t- T)At,J 



so that Gronwall's lemma finally gives the claimed estimates. D 

In order to apply the abstract framework described in Section 3, we also need the 
following higher order continuous dependence estimate. 

Lemma 6.2. Let Zqi, Zo2 € "Hi be any pair of initial data so that /xoi = /(^o«) ~ A^qj ^ 
L'^{Q), i = 1,2 and let gi,g2 € L^q^(— cxd, t; L^-y(f2)) be any pair of symbols. Then there 
exists a constant C such that, ifzi{t), i = l,2 are the solutions of (1.1)-(1.2) with initial 
data Zoi at time r and symbol gj, then the following estimate holds 

|2^^ ^^Q(A.,.,i-r)(l+||zoi|||,^ + ||zo,||^^+Uill + lmli) /'||zo||^^+ f'\g{s)\lds\ , 

where Q is a nonnegative increasing function of its arguments. 

Proof. We start by multiplying the first equation in (6.1) by 2Au = — 2PAu, getting 

d_ 

dt' 

2 ((u ■ V)ui, Au) - 2 ((u2 ■ V)u, Au) - 2 (At/^i V^i, Au) + 2 (A^2VV^2, ^u) + 2 (g, Au) 



^|Vu|2 + 2z/|Au|^ 
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The product of the third equation in (6.1) by 2A'^iIj, after an integration by parts, provides 

= - 2 (u ■ Wi, A^^) - 2 (u2 ■ V^, AV> - 2 (/'(^i) VA^, VA^') 

- 2 ([/'(^i) - f (^2)] VA^2, VA^) - 2 (/" (^1) V^A^i, VA^) 

- 2 ([/" (^1) - r (^2)]VV'2A^i, VA^) - 2 (r(V;2)VV^2A^, VA^;) 

- 4 (r (^0 V'^i V^, VAtA) - 4 an^i) - f"{^2))V'^iV^2, VA^) 

- 4 (r (^2)V2V^V^2, VA^) - 2 (r (^i)VV^ I V^il' , VA^) 

- 2 (/'"(V;i)VV^2V^ ■ (V^i + VV^2), VA^) 

- 2([r(^i) - r(v;2)]vv^2ivv^2r, vA^). 

Adding together the two equations, by Assumption (A. 5) we obtain 



(l|u|lHi,,,,(n) + |A^|^) + 2z/|Au|2,^^(^) + 2|A2^|^ 



dt 

^ - 2 ((u ■ V)ui, Au) - 2 ((u2 ■ V)u, Au) - 2 {AipWipu Au) - 2 (A7/'2VV^, Au) 



+ 2 (g, Au) - 2 (u ■ Wi, A^^;) - 2 (u2 • V^^, A^^) + 4a| VA^I^ 

- 2 ([/'(^i) - /'(^2)] VAV;2, VA^;) - 2 {f"{i,,)Vi>Ai,,, VA^j) 

- 2 ([/" (V'l) - r (V^2)]VV^2A^i, VA^) - 2 {f"{ij,)VtP2A^P, VA^) 

- 4 (/"(^AOV^T/^iV^, VAtA) - 4 ((/" (^i) - n^2))V'^iV^2, VA^> 

- 4 (r (^2) V2V'VV^2, VAV-) - 2 (/'"(^i) VV^ I VV^i I' , VA^) 

- 2 (/'"(^i)V^2VV^ ■ (VV^i + V^2), VA^) 

- 2 ([/'"(^i) - r(^2)]V^2|VV^2p, VA^;) . 

We now show that all the eighteen terms on the right hand side of the last inequality 
can be bounded by 

h{\Vn\l+\A^l^\l), 

where h is an integrable quantity. Standard computations for the Navier-Stokes equation 
lead to 

1-2 ((u ■ V)ui, Au) - 2 ((U2 ■ V)u, Au)\ 
^C\u\l/"\Au\'J"\Vui\2\Au\^ + C|u2|2^'|Vu2|2^'|Vu|2^Vu|2^' 
^-\Au\l + C|Vui|^|u|^ + C|U2|2|VU2|2|VU|^. 
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Exploiting Agmon's inequality and the interpolation inequality |0|oo ^ C*|0|2 l|0||/f-3) '^^ 
can bound the following three terms 

I -2 ( AV^ V^i , Au) - 2 ( A^2 V^, Au) + 2 (g, Au) \ 
^2[|A^|oo|VV^i|2 + lAV^ahlVV^loo + |g|2]|Au|2 
^2[\A^p\l/'\A'^p\l/'\VH2 + lA^^^M V^|f lA^I^^^ + |g|2]l^u|2 
^'^\Au\l + l^\A'^p\l + C\V4j,\t\A4;\l + C| A^^I^IV^I^ + C\g\l 

as well as the next two 

|-2 (u ■ V^i, AV> - 2 (u2 ■ V^, A^^) I 
^2|u|oo|VV^i|2|A2V;|2 + 2|u2|2|VV'|oo|A2^|2 
^C\u\l^'\Au\'J'\VH2\^'i^\2 + C|u2|2|V^|f lA^^lf 
^^|Au|^ + llA^v^l^ + C\VH>\1 + C\u2\l\Vi;\l. 

The terms arising from the double well potential can be treated using similar techniques. 
By interpolation, we have 

4a|VA^|^^i|AV|^ + C|A^|2, 

while, by (6.4), we obtain 

|-2([/'(7Ai)-/'(V^2)]VAV;2,VA^)| 

^2|/'(V'l) - /'(V^2)U|VA^2|p+3|VA^|p+3 

p+i 
^C (l + |F(V^Olf ''"'' + \F{^2)\f''^'') \V^\,\A'ij2\2\^'ij\, 

^^|AV|^ + C (1 + |F(^i)L + |F(^2)li)' |AV2|^|V^|i 

Using also Korn's inequality, from Assumptions (A. 2) and (A. 3) we deduce 

|-2 (/" (^i) V^Az^i, VAz^) - 4 (r (^i) V'^i V^, VAt/^) | 

^C|r(^l)|H±3|V^|p+3|AV^l|p+3|VA^|p+3 

p 
^C(l + |F(^Olf ''"'OlAV^MA^iUIA^^b 

^^|A^^|2 + C(l + |F(V;i)|i)^ |A^i|^^.3|A^|^. 



Arguing as in the proof of Lemma 6.1, namely, exploiting (6.4) for the first two terms, 
and Assumptions (A. 2) and (A. 3) for the second two, as well as Korn's inequality again. 
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we obtain 

1-2 ([/" (^i) - /"(t/^s)] V^aA^i, VA^) - 4 ((/" (^i) - /" (^^2)) V^i^iV^2, VAV^> 
- 2 (/" (^2) V^sAV^, VAtA) - 4 (r (V^2) V W^2, VA^;) | 
^C (1 + |F(^Olf ^^''"'^ + |F(V^2)ir'*/'^"^^) |AV^2|2|AV^iU3|VV^|2|A2^|2 
+ C (1 + |F(V;2)|f ^^'^) \Aij,\,\VAi^\,\A'i:\, 



1 



^-lA^V'l' + C (1 + |F(^i)|? + |F(^2)|?) |A^2|^|AV^i|J+3|V^| 



■7 



+ C(l + |F(^2)lt)|A^2|^|A^|i 
We are left to consider 

(6.6) I -2 (^(7/;l)VV'|V^^l|^ VA^) - 2 (r'(7/>i)V^2 V^ ■ (W^i + V^2) , VAt/;) 

- 2 ([/'"(V^i) - r (^2)]V^2|V^2p, VA^;) I 
^C(l + |F(V'i)ir^'/''' + ^^0(|AV^i|^ + lA^iy |A^2|2 + |AV;2|^)|AV;|2|AVl2 
+ C (1 + |F(^i)ir ^^Z'''^^' + |F(^2)ir''/^^"^') |AVy^|V^|2|AV|2 

^^lA^Vl^ + C(l + |F(V^i)|?)(|AV^i|^ + |AV^2|^)|AV'|^ 



C{l + \F{ij,)\l + \F{^p,)\l)\Aij,\l\Vi^\ 



Remark 6.1. In the case p E [1,2), under the assumption f^™\y) bounded for y G M we 
can still derive the estimate for the term 

(ir(^i) -r(^2)|V7A2|V7A2p,VA^) ^ C|VV^|2|A^2|^|VA^|^+3, 

which gives the same result as above. 

From the above inequalities, collecting terms we obtain 

j^{\Vvi\l+\Ai,\^^ +v\Au\l + \A'^\l 

^C\g\l + C (iV^il^ + iVuil^ + |u2|^|Vu2|^) |Vu|2 

+ C {\u2\l + \Aij2\l + (1 + |F(^i)|? + \F{i^2)\V> (IAV2I2 + |A^2|^ + |A^2|^|A^i|^^3)) \Vi^\l 
+ C(l + \W^M + (1 + |F(V;i)|?) {\A^M + |AV^2|^ + IA^il^+3) + (1 + |F(^2)lt) |A^2|^) lAV^I^. 

Denoting by 

/l = C (1 + IVV^il^ + |U2|^ + IVUil^ + |U2|^|VU2|^ + |A^2|2 

+ (1 + |F(^i)|? + |F(^2)|?) (lAVal^ + |A^2|^ + |A^2|^|A^i|J+3) 
+ (1 + WmW) {\Mi\l + |A^2|^ + |A^i|J^3) + (1 + |F(V;2)|f) |Az^2|^) , 
the above differential inequality reads as 

(6.7) ^l|z(t)||^, ^Mt)N(t)||?,,+C|g(t)|^ 
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and depends on p only through the constants C included in the definition of h. Therefore, 
the estimate obtained by Gronwall's lemma below has the structure 

where, for any potential F satisfying the assumptions (A.1)-(A.7), the function Q depends 
on At^T only through some of its powers, which, in particular, do not depend on the growth 
exponent p (i.e. the shape) of F. D 

7. Time regularity 

In this section we evaluate the distance in Tio between the solution and the initial 
datum in terms of the time span (Lemma 7.1), and we show a smoothing property for 
difference of solutions (Lemma 7.2). This will be crucial to show that Assumptions (H.l) 
and (H.3) in Theorems 3.2 and 3.3 hold true for system (1.1). 

Lemma 7.1. Given any symbol g satisfying (B.l) and (B.2), there exists a positive 
constant C such that the solution z(t), departing at time r from an arbitrary initial datum 
zq G Til, so that /io G L'^{Q) satisfies 

Mt)-Mno 

<Vt^ (llzoll?,, + l^ol^ + Al + (t - r)A,0 « + it- ^)Kr) e^(^-+(*-Mt), 

for r ^ t ^ to ■ 

Proof. The different features of the Navier-Stokes and the Cahn-Hilliard equations force to 
handle separately the two variables. We preliminarily observe that, denoting the solution 
and the initial datum as z(t) = (u(t),'?/;(t)) and Zq = {uq^i/jq), respectively. 



|u(t)-uo|2^ / \dtu{s)\2ds ^ Vt-r \\dtu\ 



L2(r,t;L2;^{n)), 



meaning that we only need to properly bound the last norm. The product of the first 
equation in (1.1) by 2dtU gives 

(7.1) '^^l^'^la + 2|<9tu|2 = -2 (u ■ V)u, dtu) + 2 (/xV^, dtu) + 2 (g, dtXi) . 

Here, having observed that 

\2{fiVij,dtu)\ = \2{tpVfi,dtu)\^2\Vi2\Moo\dtu\2^C\Vfi\MT\AtP\l^'\dtu\^, 
the right hand side can be controlled as 

1-2 ((u ■ V)u, dtu) + 2 (/iV^, dtM) + 2 (g, dtu) I 
^\dtu\l + C|u|2|Vu|2|u|H2_^^^(f,) + ClV/xl^lVblAV'ls + C\g\l. 
Replacing this estimate in the differential equality above, we have 
d 



J^^|Vu|2 + \dtu\l ^ C (|u|2|Vu|2|u|Hg^^^^(Q) + iV/il^l^hlA^b + \g\i) 
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thus, integrating in time over (r, t), thanks to Lemmas 4.2 and 5.4, we deduce 



t 

2 



\dt\i{s)\^ds 

<CA,^ (llz(r)ll?,^ + \^^{r)\l + Al + {t - r)A,) (^^ + (t - r)Al) e^(4.H-(*-)40, 

which provides the desired estimate. We now turn our attention to the order parameter. 
By interpolation, exploiting (5.11) and Lemma 5.4 again, we obtain 

||V'(t,r)-^o|lii(n) 
^|?/^(t, r) - ^o\2\\^{t, t) - ^o\\m{n) 

<:CVt^^( [ \dMs,T)\ldsY' sup ms)\\HHn) 

<,C^t^r (||z(r)||^^ + \^{r)\l + ^^ + (t - r)A,,) (A^ + (t - r)4J e^K^+(*-)^t).n 

The following smoothing property is crucial to show that our problem fits in the theo- 
retical setting of [20], which was presented in Section 3. 

Lemma 7.2. There exists a positive function Q{-, ■), increasing in both arguments, such 
that, given a pair of symbols gi, g2 satisfying (B.l) and (B.2) and any pair of initial data 
Zoi;Zo2 G "Hi so that fiQi G L'^{Q), i = 1,2, there holds 



(t-r)llz 



2 



^[IMuo+Cil+t-r) I \gis)\lds]e 



2 I r^f^ I + ^\ |„/',m2 j„ \ „Q(vit,r,i-r)(l+||zoi||2 +||zo2|l?, +|M0ili+lM02lil 



where Zj(t) stands for the solution to problem (1.1)-(1.2) corresponding to symbol gi orig- 
inating at time r from the initial datum Zjo- 

Proof. Multiplying (6.7) by (t — r), we obtain the differential inequality 

^((t-r)||z(t)||^J 
^||z(t)||^^ + C{t - r)|g(t)|^ + h{t){t - r)\\z{t)\\l^, 

where the function h is given as in the proof of Lemma 6.2. By the second estimate in 
Lemma 6.1, we deduce 

I' MsWn^ ds<:C (||zo||?,„ + £ |g(.)|^d.) (a:;" + (t - r)A;') e^(^^+(-^)^^), 
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while the integral of h can be bounded as in Lemma 6.2. Thus, the Gronwall's lemma 
entails 

(t - r)\\z{t)\\l^ 
^C (||zo||?,„ + (1 + t - r) ^' |g(.)|^ d.) {A^ + it- t)A;') e^(^^+(*-)^^ 

Q(At,T,t--r)(l+|lzoi|||^ +||zo2|||i +|M0il2 + lM02li) 

which is, the desired estimate. D 

8. Proof of the main results 

In this section we show how, properly choosing the spaces and the operators, relying on 
the results of previous sections, we can apply Theorem 3.2, and subsequently Theorem 3.3, 
to our system so to prove Theorem 2.1 and Corollary 2.2. 

Let V and H be the spaces "Hi and "Ho respectively. Observe that, whenever the 
symbol g satisfies assumptions (B.l) and (B.2), thanks to Theorem 4.1 and Lemma 6.1 
the solution operator associated to system (1.1)-(1.2) is well-defined and continuous on 
H. Moreover, thanks to Corollary 5.5, in studying the asymptotic behavior of solutions 
of (1.1)- (1.2) we can further restrict our attention to the bounded subset of V given by 

B = {zeV\ \\z\\n, + |/i|2 ^ C'g(to)}, 

which is uniformly (w.r.t. the diameter of the set of initial data) absorbing for the solutions 
of (1.1)-(1.2). Since the constraint \fi\2 = \f{'ip) — ^'ip\2 ^ C is closed w.r.t. the topology 
of V, we can further restrain our attention to the set 

when discussing the exponential decay of solution towards an exponential puUback at- 
tractor. 

Let r G M be given and let g satisfy (B.l) and (B.3) (so that (B.2) holds true as well 
for t ^ to; which is enough for our scopes), we denote by Ug{t,T) the solution operator 
to problem (1.1)-(1.2) at time t ^ t with symbol g and initial data in V. Thanks to 
Lemmata 5.4 and 6.2 the process Ug{t,T): V — )■ V^ is well-defined and continuous on 
O'^. Therefore, if to is the time appearing in assumption (B.3), the restricted family 
{f/g(t,r): T ^ t ^ to} belongs to the class W(V,to). Since the set (9^ is uniformly 
absorbing for the family {f/g(t, r)}, fixing 6=1 and the time span 

ro = 5 + T(|Of(i?)|), 

we deduce by Corollary 5.5 and Lemma 7.2 that Ug{t, t— ro) G Si^iiB), where the constant 
L depends increasingly on To, C*g(to) and Mg(to). 

In this framework we can now verify assumptions (H.1)-(H.3) for system (1.1)-(1.2). 
Indeed, (H.2) is a straightforward consequence of Lemma 6.2. The coupling between 
the Navier-Stokes and the convective Cahn-Hilliard equations makes more involved the 
validation of assumptions (H.l) and (H.3). In order to avoid further requirements over 
the symbols but (B.3), we exploit Lemma 7.1 and interpolation, thanks to a smoothing in 
the solution. The technical details of our argument are contained in the following lemma. 
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Lemma 8.1. Assume that g satisfies (B.l) and (B.3). Then there exists a positive con- 
stant C depending on the exponent q in (B.3) such that, for any initial datum zq G Oi{B), 
the solution z{t) = {u{t),tlj{t)) = ?7g(t,r)zo satisfies 



(8.1) ||u(t)||jj(.,-.,/ +||^(t)||^a^C, 

0,div ^ ■^ 

for any t ^to — 1 and r ^ t — 1 — tq . 

Proof. In this proof, we consider the two equations separately: first, as in [20], we apply 
the Giga-Sohr argument (see [15]) to the equation 

(8.2) dtu - z/PAu = h, 

where 

h = -Pu ■ Vu + P/iV^ + g. 
By Holder's and Gagliardo-Nirenberg's inequalities, we obtain 

|u- Vu|2^C|u|,|Vu|2,/(,_2) ^C|Vu|2*'"''/'|Au|2/'. 
Moreover, recalling that f{ii))Vi' G L^(i7) and f{ilj)'Vip G L^j^(l])^ 

|P/iV^|2 ^ \AiPVij\2 ^ C|AV'|^"''/''|AV|2^^ 

(actually this estimate is not optimal but, due to the previous estimates on the velocity 
field above, this does not have any influence on the final outcome) then Corollary 5.5 and 
assumption (B.3) ensure 

/ |h(s)|^ds^g(Mg,,(to)), r<:t-6-T{\0'i{B)\)=t-l-To, 

Jt~2 

which is identical to [20, Equation (84)]. From this estimate, arguing as in [20] we deduce 
the bound on u 

l|u(t)||H<=.^-)A(f,)^Q(^g,,(^o)), 

which is the first part of (8.1). 

We now turn our attention to the Cahn-Hilliard equation. The product of the third 
equation in (1.1) by —2Adtip leads to 

^|VAV'|2 + 2|VatV|^ 

= - 2 (u ■ Vdt^P, A^) + 2 (/'(^) VA^, VdtiP) + 2 (/"(^) VMV^, V^tV^) 

+ 4 (/"(V;) VV^VV, V9i^) + 2 (/'"(^) V^l VV^|^ Vdtij) . 

By Ladyzhenskaja inequality and interpolation, the first term on the right hand side can 
be controlled as 

\-2{u-VdtiP,AiP)\ 

^2|u|4|A^|4|V<9j^|2 
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and, similarly, the second one is 

|2(/'(^)VAV^,V9tV')| 
<:2\f'{i:)U\VAijU\Vdtij\, 

Taking advantage also of Agnion's and Korn's inequalities, we compute 

|2 {f"{ij)VijAij, VdtiP) + 4 (/"(V;)V^VV, V9iV^> I 
^2\f'{^)U\Vilj\oo\AijU\Vdtij\2 + Mf'm4\VijU\V^^ljU\Vdt^\2 

^c\f"m4\'^i^\T\AiJ\2\A'^\T\vd,i^\, 



12 



^\\vd,^p\i + lA^^i^ + c\f"mt\'^^\i\^^\t- 

Finally, 

|2 (r(^)v^ivv^p, vdt^)\ ^ c\rmA^^\i\^dt^\2 ^ \\vdtij\i + c\rmi\A^\i. 

Collecting the above estimates, we obtain 

^iVAi^ll^h, 
at 

where 

h = C (lA^V^I^ + |/'(V')|^|A^|^ + |ur/^|Vu|f |A^|^ + \f"ml\Vi^\l\A^\l + \rml\^^\\ 
Having observed that by Assumption (A. 6) we have 

ifwu+irwu + irwu^i/wu, 

then Corollary 5.5 provides 

/ /i(s)ds < Q(Mg(to)), r^t-5-T(|Of(fi)|), 

by the Uniform Gronwall's lemma, we deduce 

|VAV;(t)|^^g(M,(to)) t^r + 5 + r(|Of(i?)|). D 

We can now show that (H.l) holds true. Having fixed g satisfying (B.l) and (B.3), we 
denote by C a generic positive constant depending only on Mgg(to)- Then we observe 
that by Corollary 5.5 there holds 

sup ||t/g(t, r)zo||Hi ^ C 

zoGOf(B) 

for t ^ to and t ^ t — tq. Besides, Corollary 4.3 entails 

(8.3) A,r^C, t^to, r^t-To, VzoGOf(5). 
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Therefore, arguing as in Lemma 7.1, for any initial datum zq G Oi{B), we have 
(8.4) ||[/g(t, r)zo - f/g(t - s, r)zo||wo ^ C^, 

for t ^ to — 1) ^ s ^ 1 and r ^ t — tq. To proceed in our argument, we need to 
consider the two variable separately: thus, for any initial datum Zq G 0^{B), we set 
(u(t,r),'z/'(t, r)) = Ug{t,T)zQ. Provided that t ^ to — 1, ^ s ^ I, Tq ^ r ^ 2ro, by 
interpolation and Lemma 8.1, we have 

||u(t, t-r)-u{t-s,t~s-r) ||hi_^;^(o) 

^C'||u(t,t - r)-u{t-s,t-s- r)\\'l;';^;^-'\ 
as well as 

Wipit^t-r) -ip{t- s,t- s-r)\\H2(n) 
^C\m, t-r)-ijit-s,t-s- r)\\f,^^^\m, t-r)-ijit-s,t-s- r)\\'^\^^ 
^C||V^(t, t-r)-^Pit-s,t-s- r)\\'/;,^^y 
Thus we are left to control a "Ho-norm which can be split into two parts as 

||t/g(t, t - r)zo - f/g(t - s,t- s- r)zo\\-Ho 
^||[/g(t,t- r)zo - f/g(t- s,t-r)zo||-Ho + \\Ug{t - s,t - r)zo - Ug(t - s,t - s-r)zo||Ko 
^C^/s+ \\Ug{t- s,t-r)zo - Ug{t - s,t - s - r)zo\\no, 

thanks to (8.4). In order to control the last term, we observe that, by (4.7) with initial 
datum f/g(t — r, t — s — r)zo, it follows 

= (l + ||z(t -r,t-s- r)\\l^^ + 2\F{ij{t - r,t - s - r))\i + J \g{s)\l dsj 

^C ((1 + llzoll?,.^ + 2|F(^o)|i)e^^^ + Mg(to)) 

^C(l + Mg(to)). 

Exploiting Lemma 6.1 with Zi = Zq and Z2 = Ug(t — s,t — s — r)zo together with Tq ^ 
r ^ 2tq and (8.3), the desired norm can be written as 

||t/g(t- s,t-r)zo - f/g(t- s,t- s -r)zo||-Ho 
= ||f/g(t — s,t — r)zo — Ug(t — s,t — r)Ug{t — r,t — s — r)zo||-^o 

^C||zo - f/g(t -r,t-s- r)zo|k„e^(^+*^^(*°)'). 
Finally, we can use Lemma 7.1 and (8.3) so to obtain 

||zo - f/g(t -r,t-s- r)zo\\no ^ C^ {\fi{t -s-r)\l + A?_^_j_,_^) v4^_^^t_,_^e 
where the first term on the right hand side is bounded by 

Ht - r - s)\l ^ Clito) 



C^t-r,t-s-r 
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as a consequence of Corollary 5.5 and of the absorbing set considered. This together with 
the Holder's inequality and (8.3) yields the uniform estimate 

A-r,i-s-r ^ C(l + Mg(to)) 

for all t ^ to — 1, ^ s ^ 1 and Tq ^ r ^ 2tq. We thus obtain 

||zo - f/g(t -r,t- s- r)zo||^o ^ C^/s, 
which, replaced in the above inequalities, gives 

||t/g(t,t-r)zo - f/g(t- s,t- s -r)zo\\ni 
^C\\U^it,t- r)zo - f/g(t - s,t - s - r)zoC;"'^'"'-'' 
^C<.'-^V4(,-i) + c'||zo-f/g(t-r,t-s-r)zo||i;^'/^'-^^ 

for any Zq G Oi{B) and t ^ to — 1, ^ s ^ 1 = e^, Tq ^ r ^ 2tq, so that (H.l) holds true. 
We now turn our attention to (H.3), whose proof is now straightforward and follows 
from interpolation. Lemmata 8.1 and 7.1: indeed, for any zq G B, 

\\U^{t,t - r)zo - Ug{t - s,t - r)zo\\n, < Cs"'~''/''''-'\ 

for any t ^Iq, Tq ^ r ^ 2ro, ^ s ^ 1 = eo. This end the proof of Theorem 2.1. 

Finally, if (B.3) holds uniformly for to ^ I^ (exactly as (B.2)), then it is easy to see 
that the above argument applies for all times to the process Ug{t, r). In particular, (H.4) 
follows from Lemma 6.2. This proves Corollary 2.2. 
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